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We present the results of simulations of forced turbulence in a slab where the mean kinetic helicity has a maximum near the mid-plane, 
generating gradients of magnetic helicity of both large and small-scale fields. We also study systems that have poorly conducting buffer 
zones away from the midplane in order to assess the effects of boundaries. The dynamical a quenching phenomenology requires that the 
magnetic helicity in the small-scale fields approaches a nearly static, gauge independent state. To stress-test this steady state condition 
we choose a system with a uniform sign of kinetic helicity, so that the total magnetic helicity can reach a steady state value only 
through fluxes through the boundary, which are themselves suppressed by the velocity boundary conditions. Even with such a set up, the 
small-scale magnetic helicity is found to reach a steady state. In agreement with earlier work, the magnetic helicity fluxes of small-scale 
fields are found to be turbulently diffusive. By comparing results with and without halos, we show that artificial constraints on magnetic 
helicity at the boundary do not have a significant impact on the evolution of the magnetic helicity, except that "softer" (halo) boundary 
conditions give a lower energy of the saturated mean magnetic field. 
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1 Introduction 



Stars with outer convection zones tend to possess magnetic fields that display spatio-temporal order with 
variations that are often cyclic and, in the case of the Sun, antisymmetr ic with respect to the equatorial 



plane . Simulations now beg in to reproduce much of this behavior (see, e.g.. lBrown et adl2010l . lKapvla et al. 
201CL iGhizaru et a/.ll2010l ). A useful tool for understanding the outcomes of such models is mean- field 



dynamo theory. A central ingredient of this theory is the a effect. This ef fect quantifies a component of the 
mean electromotive force that is proportional to the mean mag netic field (|Moffattll 19781 . iKrause and Radlerl 
1980h . 

Mean-field theory gives meaningful predictions w hen to expect cyclic or steady solutions, and what the 
symmetry properties with respect to the equator are ( Brandenburdfl998l ) . Even in the nonlinear regime, the 
simple concept of a quenching, which reduces a locally via an algebraic function of the mean magnetic field, 
tends to give plausible results. However, under some circumstances, it becomes quite clear that this simple- 
minded approach must be wrong. Such a special case is that of a triply-periodic domain. Astrophysically 
speaking, such a model is quite unrealistic, but it is often employed in num erical simulations. It was als o 
employed as the primary tool to compute a quenching from simulations (ICattaneo and Hug heslll99fil b 
These simulations suggest that a quenching would set in once the mean field becomes comparable to a 

— 1/2 

small fraction (i? m , where R m is the magnetic Reynolds number) times the equipartition value. If this 
were true also for astrophysical bodies such as the Sun, the a effect could not be invoked for understanding 
the dynamics of the Sun's magnetic field. 

Later it became clear that there are counter examples to the simple idea that a is quenched just depending 
on the local field strength. Surprisingly, simulations later suggested that even in a triply-pe riodic domain 
a larg e-scale magnetic field can be generated that can even exceed the equipartition value (IBrandenburg 
200 ll ). However, it would take a resistive time-scale to reach these field strengths, so there was still a 
proble m. Around the same time, th e idea emerged that open boundaries might help (jBlackman and Field 
2000al lbL iKleeorin et d\ 1200(1 12002ft . This is connected wi th the fact tha t an a e ffect dynamo produces 
magnetic helicity of opposite sign at large and small scales ( Seehafer 1996 . Ji 19991 ). The magnetic helicity 
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at small scales is an unwanted by-product that can feed back adversely on the dynamo. The resistively 
slow saturation phase in periodic-box simulations can then be understood in terms of the time it takes to 
dissipate this small-scale magnetic helicity. It is indeed a particular property of triply-periodic domains 
that magnetic helicity is strictly conserved at large magnetic Reynolds numbers. A possible remedy might 
then be to consider open domains that allow magnetic helicity fluxes. 

The first simulations with open domains were not encouraging. While it was possible to reach satura- 
tion more quickly, the field w as found to level off at a value that becomes progressively smaller at la rger 
magnetic Reynolds numbers ( Brandenburg and Dobler 2001 . Brandenburg and Subramanian 2005al ) . A 
possible probl em with these simulations migh t be the absence of magnetic helicity fluxes within the do- 
main. Indeed, Brandenburg and Dobler ( 200ll ) considered a kinetic helicity distribution that was approx- 
imately uniform across the domain, so there were no gradients except in the immediate proximity of 
boundaries, where boundary conditions on the velocity prevent turbulent diffusion . The situation im- 
proved dramatically when simulations with shear were considered (|Brandenburdl2005l . iKaovla et aZll2008l . 
Hughes and Proctor 20091 ). Shear provides not only an additional induction effect for t he dynamo, but it 
migh t also lead to an additional source of magn etic helicity flux within the domain (jVishniac and Choi 
2001 . Subramanian and Brandenburgl 20041 . 120061 ) . More recently it turned out that, even without shear, 
diffusion down the gradient of small-scale magnetic helicity could, at least in principle, help avoid van- 
ishing ly small saturation levels of the mean magn etic field when the magnetic Reynolds number becomes 
large (|Brandenburg et aZ.l[2009l . iMitra et al.ll2010h . 

An important goal of the present paper it to revisit this issue using direct simulations of turbulent 
dynamos without shear, and even with the same sign of magnetic helicity everywhere, but with a spatial 
modulation of the helicity within the domain. In other words, the level of turbulence is maintained at a 
high level throughout the domain, but the amount of swirl diminishes toward the boundaries. In most 
of the simulations we include a turbulent halo outside the dynamo domain where the Ohmic resistivity 
is enhanced. This might be important as several simple boundary conditions such as pseudo-vacuum (or 
vertical field) conditions fix the value of the magnetic helicity artificially, and if fluid motions through the 
boundary are prohibited, turbulent transport there is impossible. 

Our simulations also allow us to make contact with nonlinear mean-field phenomenology where the 
evolution of the small-scale magnetic helicity is taken into account. This leads then to an evolution equation 
for an additional contribution to the a effect, «m- This approach is referred to as dynamical a quenching. 
In the present paper we will also attempt to as sess the validity of so me of the corner stones of dynamical a 
quenching. Firstly, there is the magnetic a of Pouquet et al\ (jl976l ). where the fluctuating magnetic field 
generates an otM that is proportional to the current helicity of the fluctuating field. This om counteracts 
the kinetic a, and so saturates the dynamo. Secondly there is magnetic helicity conservation which notes 
that the total magnetic helicity is nearly conserved under common conditions, and so the helicity in the 
fluctuating field can be related to the helicity in the large-scale field. Finally, there is the assumption 
that the mean current helicity of the fluctuating field is proportional to the mean magnetic helicity in the 
fluctuating field. 

As noted above, a problematic prediction of dynamical a quenching is that rapid (exponential) growth 
of mean magnetic fields will be halted below equipartition with the turbulent energy. The export of small- 
scale helicity could provide a release from this constraint but will likely occur side-by-side with export of 
the mean field. The interplay between these effects can smother the dynamo even in the presence of small- 
scale helicity transport. Treatment of large-scale helicity transport proves significantly more complicated 
than that of the small-scale helicity, but we will draw some preliminary conclusions. 

In Section[2]we discuss the dynamical a quenching phenomenology. In Section[3]we describe the numerical 
setup of the simulations whose results are analyzed in Section [U Mean-field models of the systems are 
discussed in Section [5] and we conclude in Section [6l 



2 Dynamical a quenching 



We wish to use a mean-field approach to the saturation behavior of dynamos. In what follows our averages 
will be denoted by overbars and the fluctuating terms will be denoted by lower case symbols. In the 
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simulations we will be using planar xy averaging unless noted otherwise, so the mean magnetic vector 
potential is given by 

A(z,t) = J J A(x,y,z,t)dxdy/L x L y , (1) 

A = A + a, (2) 

so the mean magnetic field is B = V x A and the mean current density is J = V x B/fio, where po is the 
vacuum permeability. In the following we adopt units in which po = 1. Throughout this paper we use the 
expressions 'mean field' and 'large-scale field' synonymously. Likewise, we refer to the 'small-scale field' as 
the 'fluctuating field'. 

We will work in the Weyl gauge (zero electrostatic potential, i.e. dA /dt = u x B — 77J), and assum e 
that there is no mean velocity. We adopt the magnetic a prescription of iKleeorin and RuzmaikTn ( 19821 ). 
As such our mean-field theoretic equations are: 

fTS 

— = Vx(5- ?? J), (3) 

S = aB- ??t J, (4) 
a = c*k + «m, (5) 



LL± ~ kf-— ~ k 2 ^ 

3 Hop f 3 fi P f ^ q 



«M = ^— ; ^ k { -— - ~ k { -^a ■ b, (6) 

where E = u x b is the mean electromotive force, «k is the kinetic a effect, B 2 q = pu 2 ms is a measure of the 
turbulent kinetic energy and r] t = Tu 2 ms /3 is the turbulent diffusivity. The parameter k{ is the wavenumber 
of the energy carrying scale of the turbulence. This q? m is taken to be the back-reaction component of a 
when it is split into kinetic and magnetic components ( Pouquet et al\ 19761 ) . 

Magnetic helicity conservation can be seen from the time evolution equation of the magnetic helicity 
density = A ■ B, 

^ = - 2r] J.B-V-F T , (7) 

where J-^ is the magnetic helicity flux. The subscript T refers to total field, which is composed of mean (m) 
and fluctuating (f) fields. In systems where the flux of magnetic helicity can be neglected (such as spatially 
homogeneous systems), and when the magnetic Reynolds number R m is large enough (and therefore rj 
small) the magnetic helicity will be nearly conserved. We define the large-scale and small-scale helicities 

as 

h m = A-B, (8) 



h{ = a ■ b = hi — h m . (9) 



Averaging Equation ([7|) we arrive at: 



dh<T dh m dhf 



dt dt dt 
dh m 



+ -^ = -2 V J-B-2 V j-b-V-F T , (10) 



o/ -2£ B-2 V J B-V (11) 
^ = - 2 £-B-2 V T r b-V-F i . (12) 
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In the spirit of mean-field theory, we will scale the fluxes to gradients of mean quantities. We consider here 
only a diffusive helicity flux of the small-scale fields, 

Ft KfVTif, (13) 

while the flux of large-scale helicity will be discussed in greater detail in Section 0J In view of Equation (|6|) , 
the evolution of /if is basically equivalent to the evolution of ccm- 



3 Numerical setup 

In this paper we present both direct numerical simulations and mean-field calculations. In both cases we 
use the Pencil CodeQ, which is a modular high-order code (sixth order in space and third-order in time) 
for solving a large range of different partial differential equations. 

We consider models with and without a halo. In both cases the horizontal extent of the domain is L x x L y 
with equal side lengths L x = L y = L, with periodic boundaries. In cases without a halo the vertical extent 
of the domain is L z = L while in cases with a halo we choose L z = 2L. In the following we measure length 
in units of fej -1 , where k\ = 2ir/L is the minimal horizontal wavenumber. We will define our magnetic 
Reynolds number as 

Rm = Urmahh, (14) 

where we assume that the wavenumber of the forcing is also the wavenumber of the turbulence, and use 
a scale separation ratio kf/k\ = 3. 

At the top and bottom, we impose stress-free velocity conditions with u z = = du x /dz = du y /d z . 
At the top and bottom we impose a "vertical field" condition, A z = dA x /dz = dA y /dz. This condition 
imposes B x = B y = and hence A ■ B = 0. As this condition on magnetic helicity may be artificial, 
and the velocity boundary condition constrains turbulent transport into the boundary, we include buffer 
"halos", such that the microscopic magnetic diffusivity is given by 

J 7/0 -7T < klZ < IT 

[t/h \hz\ > 7T, 

where t/h S> t]o- We include forced turbulence, with uniform amplitude and wavenumber k{, but a relative 
helicity a = (V x / • f)/hf of 

{COsA;iz/2 — 7T < Z < 7T , . 

, 7 ~ 16) 
\k lZ \>ir. 

For detailed about t he implementation of a forcing function with variable helicity we refer to the paper by 
Haueen et ail faOOi ). 



The above system is interesti ng from a dynamical a perspective as it contains several contrasting ele- 



ments. Unlike iMitra et al\ (|2010l ). the forcing helicity is all of one sign, and so we expect the magnetic and 



current helicities to also be of one sign. This implies that the production term in Equation ([7]) be finite 
even after volume averaging, and there may never be a final steady state for the magnetic helicity. Further, 
fluxes through the actual boundaries are reduced through the velocity boundary conditions as well as the 
resistive destruction of the field in the halos. It is not clear whether the magnetic vector potential will even 
have a final steady solution. An e xample of an unste a dy m agnetic helicity in an otherwise fully steady 
dynamo was presented in Fig. 2 of Brandenburg et al. ( 20021 ). 



http : / /pencil-code . googlecode . com 
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Figure 1. Visualizations of each of the three components of B on the periphery of the computational domain, for a run with i? m = 1300 
using 256 X 256 X 512 mesh points. The lower planes show the field at the lower boundary, where it is quite weak. Note the presence of 
mild overshoot into the upper and lower halos, where rm = 250r)o- In the bulk of the domain the x and y components show a large-scale 
field with variation in the z direction, while B z does not show a systematic mean field. 



The main difference compared with earlier work is that in Mitra et al. ( 2010l ) there was an equator at 
z = with kinetic helicity of opposite sign for z < 0. Consequently, also hf changes sign, allowing an 
efficient exchange of magnetic helicity by the turbulence. The present model is more similar to that of 
Brandenburg and Doblerl (|200ll ). except that there the kinetic helicity profile was flat in the bulk of the 
dynamo interior and it dropped to zero only immediately at the boundary of the domain, or gradually so 
in those cases where a conducting non-turbulent halo was included. 

The helicity in the halos, as noted above, will be suppressed by the low conductivity, as we expect 
small-scale helicity transport away from the active central region. A strong, rapidly achieved final mean- 
field would indicate that flux of small scale helicity provides a clear escape from dynamical a quenching. 
Finally, a clear difference between halo simulations and simulations without a halo will be evidence that 
the boundary conditions are generating artificial constraints. If these differences are visible in the field 
itself (as opposed to the vector potential), they are likely due to the reduced turbulent diffusion into the 
boundary. 

Mean quantities are calculated from time series over a long stretch of time where the relevant quantities 
are approximately stationary in the statistical sense. We use the time series further to calculate lower 
bounds on the error bars as the maximum departure between these averages and the averages obtained 
from any of the three equally long subsections of the full time series. 



4 Results of simulations 



In the statistically steady state, the magnetic field shows a large-scale magnetic field that varies in the z 
direction; see Figured! It is therefore meaningful to describe the dynamics of this large-scale field by using 
horizontal averages as noted in Section [2J We will use angular brackets and capitals to refer to volume 
averages over the volume V = L x L y L' z , where L' z = Z2 — z\ and —z\ = Z2 = 2k^ 1 . This will mark the 
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Figure 2. Magnetic helicity density of the small-scale magnetic field in the Weyl gauge, h t , as a function of z and t (left-hand panel), 
and the gauge-independent magnetic helicity of the large-scale field, H m in units of 2M m /fci (right-hand panel), for Run H3 with 
7/h = 100r;o- For comparison, the Weyl-gauged magnetic helicity of the large-scale field is shown in the inset. The data are averaged over 
regular time intervals of about 0.3 diffusive times, which also explains the absence of data at t = 0. 



boundaries of a smaller domain well within the dynamo region. Thus, we write 



H(t) = {h) v V 



L x LyL' z 




h(x, y, z, t) dx dy dz = — / hi(z,t)dz 



(17) 



where h and H could stand for H m and h m , or for H{ and /if, for example. We note, however, that 
these quantities may be gauge-dependent. We also define the magnetic energy of the mean field as M m = 
(B 2 /2)V. 



4.1 Small-scale helicity flux 

We define the magnetic helicity densities for the mean and fluctuating fields as 



h 



■w 



A B, 



h 



■w 



a ■ b. 



(18) 



The superscript W indicates that we are working in the Weyl gauge modulo possible influences of the 

— w 

boundary conditions that have been mitigated through the use of halos. It turns out that h m has a 
systematic variation in time while h f does not; see Figure [2J It makes therefore sense to average the 



evolution equation for /i f in time, so we have (iMitra et a/Jl2010h 



dh 



w • 



dt 



= -2(5 • B) T - 2n(j ■ b) T 



'-^f )T, 



(19) 



where subscripts indicate time averaging over the interval T. In the Weyl gauge the magnetic helicity flux 
of the small-scale field is given by IPY = e x a, where e = E — E is the electric field for the fluctuating 
quantities. Given that the first two terms on the rhs of Equation (|19p are gauge-invariant, (V • J-{) must 
also be gauge-invariant, so we can drop the superscript W and note that, in the particular case at hand, 
we have (V • ^£ V )t = (V • J^{)t- We emphasize that (V • T{)t is still a function of z. 
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4.2 Large-scale helicity flux 

In order to assess the full magnetic helicity budget, we also need to take the magnetic helicity of the mean 

— w — w 
field into account. Since h m is time-dependent, it is not possible to invoke a similar argument as for /i f . We 

are therefore forced to abandon a detailed analysis of the z dependence of the magnetic helicity budget and 

restrict ourselves to the analysis of the volume-integrated magnetic helicity, H m , and its co rresponding flux 

divergence, Q m , using the gauge- invariant prescription of Brandenburg and Dobler ( 200ll ) for the volume 

in zi < z < Z2, 

H m (t)= f 2 hZ{z,t)dz + A( Zl ,t)xA(z 2 ,t), (20) 

and 

"22 



Q m (t) = -[E(zi,t) + E(z 2 ,t)) ■ [ 2 B(z,t)dz, (21) 



where E = r]J — £ is the mean electric field expressed in terms of horizontally averaged Ohm's law. Note 
that H m and Q m obey the evolution equation 



dHn 

dt 



2 2 Z-22 

S-Bdz-2r] I J-Bdz-Q m . (22) 

Zl J 2i 



— w . _ 

It turns out that, unlike h m and its volume-integral, H m is statistically steady (see Figure [2D, so we may 



now also average Equation (|22p over time. 



4.3 Magnetic helicity budgets 

In Table [H we summarize the helicity budgets, namely the six terms on the rhs of Equations ()19p and (|22p . 
of which the 2(£ ■ B) term occurs twice. We have used here the more descriptive symbol (V • ^F m )vT = 
Qm/( z 2 — z i) for the flux divergence of the helicity of the mean field. In order to simplify the notation, we 
drop from now on the subscripts VT and define angular brackets without subscripts as combined averages 
over a long enough time span and over the volume V in z\ < z < Z2, where again, —z\ = z 2 = 2k^ 1 . 

Tabic 1. Summary of the volume and time averaged terms on the rhs of Equations < U9t and d22l >. normalized by T?to^e q > while (-B 2 ) is normalized 
by Sg q . Runs HI— 6 refer to systems with a poorly conducting halo with wh/vo ^ ^m> while in Runs VFl— 6 refer to systems without a halo 
and a vertical field boundary condition at \k\z\ — tt. The data for Run H6 arc given for completeness, but it has not run long enough to have 
satisfactory statistics. 



Run 




(W) 


2{£-B) 


2i](J -B) 


2 V (j ■ b) 


(V • T m ) 




HI 


20 


0.56 


-0.423 ± 0.003 


-0.068 ± 0.000 


0.408 ± 0.002 


-0.360 ± 0.006 


0.018 ±0.013 


H2 


50 


0.33 


-0.208 ±0.003 


-0.018 ±0.000 


0.190 ±0.001 


-0.192 ±0.005 


0.012 ±0.004 


H3 


140 


0.15 


-0.086 ± 0.003 


-0.003 ± 0.000 


0.078 ± 0.001 


-0.079 ± 0.005 


0.005 ± 0.001 


H4 


270 


0.12 


-0.047 ± 0.002 


-0.001 ± 0.000 


0.041 ± 0.000 


-0.046 ± 0.001 


0.003 ± 0.000 


H5 


520 


0.08 


-0.024 ± 0.000 


-0.000 ± 0.000 


0.020 ± 0.000 


-0.024 ±0.001 


0.002 ± 0.000 


H6 


1280 


0.08 


-0.029 ± 0.023 


-0.000 ± 0.000 


0.009 ± 0.000 


-0.007 ± 0.007 


-0.007 ± 0.004 


VFl 


10 


0.62 


-0.823 ±0.011 


-0.163 ±0.002 


0.822 ± 0.005 


-0.669 ± 0.005 


-0.000 ± 0.012 


VF2 


20 


0.43 


-0.434 ± 0.004 


-0.051 ± 0.002 


0.436 ± 0.003 


-0.400 ± 0.005 


0.002 ± 0.026 


VF3 


50 


0.32 


-0.250 ±0.013 


-0.019 ±0.001 


0.247 ± 0.002 


-0.224 ±0.009 


0.006 ± 0.012 


VF4 


120 


0.28 


-0.138 ±0.009 


-0.007 ±0.000 


0.134 ±0.001 


-0.143 ±0.004 


0.001 ± 0.008 


VF5 


220 


0.25 


-0.091 ±0.002 


-0.003 ± 0.000 


0.082 ± 0.001 


-0.086 ± 0.002 


0.004 ± 0.003 


VF6 


400 


0.15 


-0.053 ± 0.002 


-0.001 ± 0.000 


0.046 ± 0.000 


-0.052 ±0.001 


0.005 ± 0.002 
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Figure 3. Time-averaged terms on the right-hand side of Equation I I19II for Run H4, {£ ■ B)x and r)(j ■ 6)t (left panel), the difference 
between these terms compared with the magnetic helicity flux divergence of small-scale fields (V ■ ^j V )t (middle panel), and the flux 
itself compared with the Fickian diffusion ansatz (right-hand panel). The fluxes are in given in units of r) t oB^ q and the flux divergence 
is given in units of kirj t oB^ q . 

The results given in Table Q] show that 2(£ • B) is balanced essentially by (V • ^F m ), because rj(J ■ B) is 
small. On the other hand, for the magnetic Reynolds numbers considered here (R m 500), the 2r](j ■ b) 
term is still quite large, and contributes mainly to balancing —2(£ ■ B) in the magnetic helicity balance 
for the flu ctuating field. The other (smaller) contribution comes from (V • This result is quite similar 
to that of Mitra et al. ( 20ld ) for the case of a linearly varying kinetic helicity profile, where it was found 
that, even though most of (£ ■ B) is still balanced by n(j ■ b), both (£ ■ B) and (V • Ti) vary little with 
R m and must eventually dominate over rj (j ■ b) as ry decreases with increasing R m . This was estimated to 
happen at R m = 10 3 ...10 4 . In the model presented here, this is not so obvious, because (£ ■ B) shows still 
a rapid decline with increasing R m . This may be a consequence of the fact that also B 2 declines still quite 
rapidly with increasing R m , which indicates that the quenching is i? m -dependent, at least for the values 
of R m considered here. 

Note that the final field strength for systems without halos tends to be higher than for systems with 
halos: turbulent transport of the mean field out of the active region plays an important role. This implies 
that the turbulent flux of magnetic helicity from small-scale fields has a weaker effect on the final strength 
of the mean field than the turbulent flux of the mean field itself. 

Note also that the total helicities H m and Hf show little difference in the two setups, suggesting that 
the artificially imposed h = constraint on the boundary is not generating spurious results. 



4.4 Magnetic helicity fluxes 

In Table [2] we collect results for the magnetic helicity flux divergence. The profile of the flux of magnetic 
helicity from the small-scale magnetic field, J-f, is reasonable well described by a Fickian diffusion ansatz. 
In Figure [3] we show the profiles of (£ ■ B) and rj(J ■ B), compare the residual 2{£ ■ B) — 2ij(J ■ B) with 
the divergence of the magnetic helicity flux, and finally compare the flux J-{ = e x a with that obtained 
from the diffusion approximation, — KfV/if. 

There are several additional points to be noted about the simulation results. Firstly, based on earlier 
results for triply-periodic domains one expects that H m and Hf have the opposite sign, which is indeed 
always the case. Furthermore, we expect that the current helicity of the mean fields, C m = (J -B), and H m 
have the same sign and that C m /H m ~ kf. This is indeed the case, except that the simulations give about 
half or less than the naively expected value for C m /H m . This indicates that the large-scale magnetic field 
is not fully helical, a potential reason for the modest mean-field saturation strength even in the presence 
of a magnetic helicity flux of small-scale fields. Likewise, one expects that Cf and Hf have again the same 
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sign. Again, this is borne out by the simulations, but the ratio Cf/kfHf is typically 3-5 times larger than 
the expected value of unity. This may well be a consequence of the presence of a finite flux divergence of 
magnetic helicity of small-scale fields. 

Finally, we find that the sign of the flux divergence of magnetic helicity density of fluctuating and mean 
magnetic fields has the same sign as the respective magnetic helicities themselves. This is generally the 
case and is well motivated by the Fickian diffusion ansatz. Given that that we find Kf ~ 0.3ry t o, we should 
expect that Qf/ffto^i-f^f is also about 0.3, but the real value is only 0.1. On the other hand QmA/to^i-^m 
varies between 0.2 and 0.7, but tends to decrease with increasing values of R m , although it remains above 
Qi/rjtokfHf and may be approaching it from above. As long as the transport of large-scale helicity has 
a larger transport coefficient than that of the small-scale helicity, we expect that the small-scale helicity 
transport will not result in larger helical mean field stren gths even though it allows for a stronger post- 
kinematic a effect ( Brandenburg and Subramanian 2005al ) 



The complication that rj(J ■ B) is not small even for the largest R m should not be forgotten. The total 
helicity being forced, and eventually even the halo 'buffer' zones will transmit information to the boundary. 
It is therefore not clear how well a Fickian diffusion ansatz is justified for the helicity of the mean magnetic 
field. 



5 Connection with mean-field models 



In order to perform mean- field simulations, we 
coefficients. A robust tool for extracting these 



need to include all the relevant turbulent transport 
coefficients from simulation is the test-field method 



of ISchrinner et al\ (|2005l . 120071 ); for applications to time-dependence turbulence see iBrandenburg et al. 
( 2008a bh . We apply this technique both to the kinematic and to the nonlinear stage using the 
so-called quasi-kinematic tes t -field method Brandenburg et all ( 2008c! ) ; for a justification of it see 
Rheinhardt and Brandenburg! ( 20ld ). In Figure 2] we show not only the values of a and rj t as determined 
by the test-field method, but also the 7 and 5 effects in the more general expression 



£ = aB + x B — rj t J + 5z x J. 



(23) 



As expected, the latter are negligible. Interestingly, we see evidence of quenching of r/t in the active region, 
even though the mean-field is well below equipartition ({B ) = 0.1B^ q ). Further, we see an approximately 
2- fold reduction both in a and in r] t . 

In Figure [5] we compare the evolution of (B 2 )/B^ q for the direct simulation with the solution of the 
corresponding mean-field model. There is excellent agreement in the final saturation level, and in both 

Tabic 2. Normalized values of magnetic helicity, current helicity, and magnetic helicity flux divergence both for small-scale and large-scale 
magnetic fields. 



Run 




hH m /2M m 


C m / kfH m 


Qm/vtakiH m 


k 1 H l /2M m 




Qf/jfto&iflf 


HI 


20 


-0.94 ±0.02 


0.54 ±0.02 


0.69 ±0.01 


0.15 ±0.00 


2.05 ±0.04 


0.21 ±0.15 


H2 


50 


-0.89 ±0.01 


0.56 ±0.01 


0.68 ±0.02 


0.21 ±0.01 


2.50 ±0.12 


0.18 ±0.07 


H3 


140 


-0.93 ±0.06 


0.51 ±0.01 


0.55 ±0.03 


0.34 ±0.00 


3.02 ±0.06 


0.08 ± 0.02 


H4 


270 


-0.97 ±0.02 


0.50 ±0.01 


0.41 ±0.01 


0.38 ±0.00 


4.29 ±0.07 


0.08 ±0.01 


H5 


520 


-0.90 ± 0.02 


0.53 ±0.01 


0.34 ±0.01 


0.46 ±0.00 


4.96 ±0.09 


0.07 ±0.01 


H6 


1280 


-1.30 ±0.01 


0.36 ±0.02 


0.08 ±0.07 


0.36 ±0.04 


6.68 ±1.40 


-0.15 ±0.18 


VF1 


10 


-2.47 ±0.45 


0.20 ±0.03 


0.45 ±0.08 


0.11 ±0.02 


2.67 ±0.52 


-0.21 ±0.36 


VF2 


20 


-2.40 ± 0.43 


0.21 ±0.03 


0.40 ± 0.08 


0.16 ±0.02 


2.69 ±0.31 


0.02 ± 0.40 


VF3 


50 


-2.26 ±0.37 


0.22 ±0.03 


0.31 ±0.07 


0.24 ±0.03 


2.68 ±0.26 


0.06 ±0.14 


VF4 


120 


-1.82 ±0.28 


0.27 ±0.04 


0.29 ±0.04 


0.28 ±0.01 


3.49 ±0.07 


0.02 ±0.10 


VF5 


220 


-1.84 ±0.32 


0.26 ±0.05 


0.19 ±0.04 


0.32 ±0.00 


3.83 ±0.07 


0.06 ± 0.04 


VF6 


400 


-1.49 ±0.23 


0.32 ± 0.04 


0.23 ±0.04 


0.43 ±0.01 


4.76 ±0.06 


0.07 ± 0.04 
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Figure 4. Profiles of a and rjt for Run H4 for the saturated case (solid lines) and the kinematic case (dashed lines), obtained using the 
test-field method. In the left and right panels, we also show 7 and <5, respectively. 
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Figure 5. Comparison between the saturation behavior in the simulation (left) and the mean-field model (right) for the same parameters: 
K m /3 = 100, C a = 3. Magnetic helicity flux and its fit in terms of a Fickian diffusion law, for Run H4. 

cases the amplitude overshoots slightly before settling at a somewhat lower value, but the kinematic 
growth rate is much faster in the mean-field model than in the simulation. This discrepancy is not yet well 
understood and should be reconsidered in future work. Perhaps significantly, the rise time of the mean field 
i s rapid even in terms of the turb ulent turnover time. This means that memory effects become important 
( Hubbard and Brandenburg 20091 ). and that the actual growth rate would be reduced compared with that 



obtained from simple estimates. The overshoot may simply be an artifact of the finite time it takes the 
system to convert mean fields into correlated small-scale fields, which is not included in the mean-field 
model. 



6 Conclusions 



Confirming earlier work of lMitra et al 



we have demonstrated the existence of a diffusive flux J-i 
of mean magnetic helicity of the small-scale field. In the present case, however, the Weyl-gauged magnetic 
helicity of the large-scale field never reaches a steady state. Nevertheless, the magnetic helicity density of 
the small-scale magnetic field is found to be statistically steady, so the corresponding magnetic helicity flux 
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must be gauge- independent ( Mitra et al. 2O10l ). This supports the validity of using the small-scale magnetic 
helicity as a meaningful proxy for the small-scale current helicity, and hence the magnetic correction to 
the a effect. 

Understanding the transport of magnetic helicity of the large-scale field, would be useful for creating 
analytic post-kinematic models. Although we have not converged on a formula for this flux, it is certainly 
finite and apparently R m dependent. It is not yet clear whether this flux will converge to a diffusive one 
for large i? m . Our mean-field simulations reproduce the final field strength well, reinforcing the conclusion 
that post-kinematic dynamical a quenching can be used as part of a mean-field simulation. 

The preliminary evidence on the use of small-scale helicity fluxes to escape the small predicted post- 
kinematic mean fields is negative: the observed flux of large-scale helicity, while poorly modeled, is larger 
than the flux of the small-scale helicity. If this holds for larger R m , it would have the unfortunate result of 
closing escape holes from a quenching opened by JFf , but would also imply that dynamo systems with more 
realistic profiles than simple homogeneity will reach R m independent behavior for high but currently nearly 
numerically achievable R m . It is likely that conclusive e vidence for or against i? m -dependent quen ching 
requires values of R m in the range between 10 3 and 10 4 dBrandenburg et odl200gLlMitra et ofj|2010l ). 
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